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ABSTRACT 

We present the results of a 135-arcmin^ search for high-redshift galaxies lensed by 29 clusters 
from the MAssive Cluster and extended MAssive Cluster Surveys (MACS and eMACS). We 
use relatively shallow images obtained with the Hubble Space Telescope in four passbands, 
namely, F606W, F814W, FI lOW, and F140W. We identify 130F814W dropouts as candidates 
for galaxies at z > 6. In order to ht the available broad-band photometry to galaxy spectral en¬ 
ergy distribution (SED) templates, we develop a prior for the level of dust extinction at various 
redshifts. We also investigate the systematic biases incurred by the use of SED-bt software. 
The fits we obtain yield an estimate of 20 Lyman-break galaxies with photometric redshifts 
from z ^ 7 to 9. In addition, our survey has identihed over 100 candidates with a signihcant 
probability of being lower-redshift (z ~ 2) interlopers. We conclude that even as few as four 
broad-band biters - when combined with btting the SEDs - are capable of isolating promis¬ 
ing objects. Such surveys thus allow one both to probe the bright end (M1500 ^ —19) of the 
high-redshift UV luminosity function and to identify candidate massive evolved galaxies at 
lower redshifts. 

Key words: galaxies; high-redshift - galaxies: statistics. 


1 INTRODUCTION 

Reionization and galaxy formation are key events in cosmic his¬ 
tory; the former depends upon, and the latter is responsible for, 
the characteristics of high-redshift galaxies. Therefore the study of 
such galaxies is a major component of the ongoing intensive inves¬ 
tigation into the early epochs of the universe. 

The most credible redshift determinations arise from spec¬ 
troscopic analysis. Although the catalog of spectroscopically con¬ 


firmed high-redshift galaxies confines to expand (e.g.,|Richard et al.| 


201 1| Vanzella et al.|201 l|[Bradac et al. 

20 12[ [Bradley et al. 

2012[ 

Ono et al.|2012l[Finkelstein et al.|2013 

[Oesch et al.|20151 

Zitrin[ 


et al.|[20f5|>, the time required to obtain spectra of faint objects 


limits the scope of spectroscopic surveys. Hence the most fruit¬ 
ful method for expanding the catalog of high-redshift galaxies is 
the dropout technique ( [Steidel et al.|1995| |1996^ , which relies on 
multiple-passband photometry. Breaks in an object’s spectrum - in 
particular, both the Lyman break (at restframe 912 A) and the 4000- 
A break - can cause it to ‘drop out’ of passbands blueward of the 
break due to absorption of its radiation by neutral gas. Comparison 
of the observed and rest-frame wavelengths of the break immedi¬ 
ately yields a crude (‘photometric’) redshift estimate. One can sub¬ 
sequently improve this estimate by fitting galaxy template spectral 
energy distributions (SEDs) to the observed multi-band photome¬ 
try. 

Multiple researchers have applied this technique to identify 
both intermediate-redshift jSteidel et al.|1999||Ellis et al.|2001|[Gi^ 


avalisco et al.|2004l|0uchi et al.|2004[rStark et al.|2009[[0esch et al. 


3 ouwens 


et al. 
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2011 

Lorenzoni et al.|2013l [Oesch et aL|2013[ 

Schenker et al.[20 

l3l|El 

is et al.|2013|(. The widest of these surveys 


To push this technique to fainter magnitudes, other studies 


have combined it with the power of gravitational lensing dEllis 

et al.[[2001 

[Richard et al.|2006[ [2008[ [Bradley et al.|2014 

[Atek 

et al.|2014| [Zheng et al.|2014[l; it was thus that Coe et al. 

20131 


identified a galaxy with photometric redshift z ~ 11 (see also 
|Pirzkal et al.|20l5| l. Lensed surveys tend to cover a smaller solid 
angle than field surveys because of their dependence on high- 
mass foreground galaxy clusters. One of the most extensive such 
projects is the Cluster Lensing And Supernova survey with Hubble 
(CLASH— [Postman et al.|20T^ , which imaged 25 clusters in 16 
filters, with integration times in each filter ranging from 1975 to 
4920 s. (See for instance [Bradley et al.||2014l jZitrin et ^|2013[ 
[Bouwens et al.||2014| ) Most lensed surveys apply more time to 
smaller solid angles; for instance, the Hubble Erontier Field^pro- 
gram devotes 140 orbits to each of six massive clusters. Still ongo¬ 
ing, this deep-imaging program has already detected a substantial 
number of high-redshift galaxies (see Atek et al.|2014| [2015[|Zhen^ 


jet al.|20l4| [McLeod et al.|2015[[Ishigaki et al.|2015[. 


^ http://WWW.stsci.edu/hst/campaigns/ 
frontier-fields 
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Table 1. Exposure times 




.. Integration Times . 


Eil- 

Clus- 


E606W 

E814W 

EllOW 

E140W 

ters 

ters 

This work 

1200 s 

1440 s 

706 s 

706 s 

4 

29 

CLASH“ 

1975 s 

4103 s 

2415 s 

2342 s 

16 

25 


“ CLASH integration times are averages from table 5 of [Postman et al.| 
j2012) . Actual exposure times vary from cluster to cluster, depending on 
the quality of previous observations. 


Despite these successes, photometric redshifts always contain 
an element of uncertainty, given the possibility of low-redshift ob¬ 
jects mimicking the colours of high-redshift galaxies ^Mobash^ 
jet al.|[2005| [Schaerer et al.|[2b07[ [Schenker et al.[[2012^ . These 
low-redshift interlopers may be red stars or galaxies with high 
equivalent-width emission lines ( jAtek et al.|201 1| >. Thus in search¬ 
ing for high-redshift dropouts, it is important to analyse passbands 
blueward of the dropout band in order to detect the flux enhance¬ 
ment due to strong emission lines (like Ha and [OIII]) which would 
indicate a lower redshift. 

To summarize, one can increase the yield of photometric red- 
shift surveys by employment of larger sample sizes, by utilization 
of gravitational lensing from massive galaxy clusters, and by a ju¬ 
dicious choice of passbands. The most massive clusters known to 
date at z > 0.3 are those identified by the MAssive Cluster Survey 
(MACS) jEbeling et al.|2001||2007|[20T0i|Mann & Ebeling|20l2| >, 
which systematically catalogued the most X-ray luminous - and 
hence the most massive - galaxy clusters. In this work we employ 
the lensing power of 29 such clusters not studied by CLASH, thus 
conducting one of the broadest lensed dropout searches to date. 
Table [T] compares this survey to CLASH; despite the longer inte¬ 
gration times and the greater number of passbands in CLASH, the 
similarity in solid angle coverage bodes well for the identification 
high-redshift candidates by our survey. 

Throughout this paper we assume a standard concordance cos¬ 
mology with Ho = 70 km s“^ Mpc”'^, = 0.3, and Ha = 0.7. 

We express all magnitudes in the AB system (jOke & Gunn|198^. 


2 DATA 

This survey analyses Hubble Space Telescope (HST) images of 
the 28 MACS clusters and one eMACS (extended MACS) clus¬ 
ter jEbeling et al.|20f^ listed in Table This sample comprises 
all MACS clusters that were not part of CLASH and for which 
HST images exist in all of the following passbands: E606W and 
E814W on the Advanced Camera for Surveys (ACS); and EllOW 
and F140W on the Wide Field Camera 3 (WFC3). All images were 
obtained through HST SNAP-shot program^ The exposure times 
(identical within each passband across all clusters) appear in Ta- 
ble^ 

The WFC3 field of view (4.65 arcmin^), being smaller than 
that of the ACS, determines the 135-arcmin^ solid angle of this 
survey. 


Table 2. Clusters surveyed, with Milky Way hydrogen column density (in 
units of cm“^) for estimating extinction 


Cluster 

uh 

Cluster 

Uh 

eMACSJ1057.5+5759 

0.56 

MACS11354.64-7715 

2.86 

MACSI0140.0-0555 

2.85 

MACSJ1621.44-3810 

1.12 

MACSI0152.5-2852 

1.51 

MACSJ1652.34-5534 

2.33 

MACSI0257.6-2209 

2.18 

MACSI 1731.64-2252 

6.29 

MACSI0451.9-t0006 

7.23 

MACS11738.14-6006 

3.74 

MACSI0712.3-t5931 

5.43 

MACSJ1752.04-4440 

3.06 

MACSI0916.1-0023 

3.25 

MACSJ2050.74-0123 

7.50 

MACSI0947.2-t7623 

2.22 

MACS12051.14-0215 

8.14 

MACSIl 115.24-5320 

0.90 

MACS12135.2-0102 

4.27 

MACSIl 124.54-4351 

2.04 

SMACSJ0234.7-5831 

3.64 

MACSI 1133.24-5008 

1.44 

SMACSJ0549.3-6205 

4.50 

MACSIl 142.44-5831 

1.77 

SMACSJ0600.2-4353 

6.24 

MACSI1226.84-2153C 

1.87 

SMACSJ2031.8-4036 

3.91 

MACSI 1236.94-6311 

1.68 

SMACS12131.1-4019 

3.00 

MACSI1319.94-7003 

1.47 




3 ANALYSIS 

One result of the opportunistic nature of the SNAP-shot program is 
that the four images of each cluster are taken at random times dic¬ 
tated by scheduling requirements. As a result the images of a given 
cluster are not, in general, aligned. In addition, the ACS plate scale 
is significantly smaller than that of WFC3. In order to facilitate 
comparison between images in different passbands, we redrizzled 
all images (using DrizzlePac^l to the pixel scale and reference 
frame defined by the F140W image for the relevant cluster. 

The presence of low-sensitivity regions (‘blobs’ - see |Dress^ 
|2014| l can complicate the analysis of WFC3 images. These artefacts 
are the result of differential reflectivity of the Channel Select Mech¬ 
anism Mirror. In several instances (discussed in Section one of 
these blobs (or another defect) in the El lOW image coincides with 
both a detection in the F140W channel and dropout behaviour in 
F814W and F606W (see for instance the final row of Fig.|^. In ad¬ 
dition, ambient light had contaminated the majority of the El lOW 
image for cluster MACSJ0712.3-1-5931 (see for instance the fourth 
and fifth rows of Fig.|^. In both of these situations, only three pass- 
bands of usable data are available. 

We then stacked the F140W and FllOW images and ran 
SExtractor ( jBertin & Amouts|1996| ) in dual-image mode, us¬ 
ing the stacked image for detection. We employed a 12-pixel rect¬ 
angular annulus for background determination and thus obtained 
a catalog of objects with isophotal magnitudes in each passband. 
Our inital catalog of SExtractor results comprise 37,809 records, 
the flux errors of which we corrected for correlated noise according 
to the prescription of |Casertano et al.] ( |2000^ . Since the WFC3 im¬ 
ages determine our segmentation, and since the WFC3 point-spread 
function (PSF) is significantly wider than that of the ACS, we do 
not perform any PSF-matching. We want to capture as much flux 
as possible from the ACS images in order to insure that the objects 
we consider are truly dropping out in the ACS bands. By retaining 
the tighter ACS PSF we allow for more accurate detection of this 
dropout behavior. 

The next task is to discriminate between galaxies, stars, and 
artefacts. To do so, we consider both the objects’ morphology 
and their colours; we describe the morphological criterion here 


^ GO-10491, GO-10875, GO-12166, GO-12884: PIH.Ebeling ® http://drizzlepac.stsci.edu 
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Figure 1. We use both morphology and spectral energy distribution to dis¬ 
criminate between stars and galaxies; these plots show our morphological 
criteria. In both panels the line corresponding to point sources (stars) is ap¬ 
parent, running diagonally in the top panel and (roughly) horizontally in 
the bottom panel. Any object below this line is more compact than a point 
source and thus must he an artefact. Hence we discard anything in the un¬ 
shaded portion of either plot. Since the star line and the galaxy cloud begin 
to interpenetrate at higher magnitudes and since HST may not resolve com¬ 
pact faint high-redshift galaxies, we place this cutoff la above the star line 
at low magnitudes and 2a below the star line at high magnitudes. Any ob¬ 
ject in the light green (upper) region of both plots (3a above the star line) we 
accept as a galaxy without further examination. The remaining objects (in 
the middle, hlue-shaded region) we subject to the spectral criterion outlined 
in Section|4j2] 


and the spectral criterion in Section |4.2[ As Fig. shows, point 
sources (stars) occupy a well-defined region (‘star line’) in both 
magnitude/surface brightness space and magnitude/half-light ra¬ 
dius space. Any object which lies more than 3a above this line in 
both plots (i.e., in the light green regions of both panels of Fig.[^ 
we accept as a galaxy without further examination. 

Ideally we could now reject anything on this line as a star and 
anything below this line as an artefact. However, at higher magni¬ 
tudes the galaxies begin to bleed into the star line; in addition, high- 
redshift galaxies can remain unresolved by WFC3 ( [Oesch et al.| 
|2010| >. To account for these facts we place the rejection limit Icr 
above the star line at low magnitudes but 2 a below the star line at 
high magnitudes; anything that lies below this limit in either plot 
(i.e., in the unshaded region of either panel of Fig. [TJ we reject as 
being either a star or an artefact. 


The remaining objects lie in the blue region of the figure; we 
provisionally admit these objects into the next phase of our analysis 
but will use the spectral criterion of Section [43] to eliminate M-stars 
and brown dwarfs. 

At this point we also exclude detections in the noisy portions 
of the WFC3 field of view|f] 

We identify F814W dropouts by requiring at least a 5a de¬ 
tection in either WFC3 band (F140W or FllOW) and less than 
2a signal-to-noise ratio in both ACS bands (F814W and F606W). 
We visually inspected the dropouts, eliminating diffraction spikes 
and areas in which a nearby bright source had corrupted the SEx- 
tractor results; we also eliminated ‘detections’ that appeared to be 
serendipitously grouped noise. 

In addition, we checked the SExtractor segmentation map for 
these sources to ensure that SExtractor properly discriminated be¬ 
tween the objects themselves and neighboring sources. For cases 
in which it did not, we set an appropriate aperture for each object 
and reperformed the photometry in those apertures. In some cases 
the aperture photometry resulted in at least a 2a detection in the 
dropout band, causing us to eliminate these sources from consider¬ 
ation. Our final / 8 i 4 -dropout catalog consists of 130 sources. 

Finally, to account for foreground (Milky Way) extinction, we 
first convert the column densities tih from Tabl ej^to values of Ay 

j2009 1 and then apply the 


Giiver & Ozel 


using the prescription of 
ICardelli et al.|(T989|> law to calculate the extinction in each band. 


4 SED FITTING 

4.1 Obtaining Redshift Probability Distributions 

Both the Lyman break and the 4000-A break can cause dropout be¬ 
haviour. Thus, a decrease in F814W flux accompanied by a strong 
detection in FllOW could reflect the 4000-A break redshifted to 
1-5 ^ z < 2 oi the Lyman break redshifted to 7 < 2 < 10. The 
lower-redshift objects are likely to be massive, passively-evolving 
galaxies. These galaxies (at such redshifts) typically are quiescent, 
extremely compact, and already quite old, with mass densities at 
least an order of magnitude greater than those of local elliptical 
galaxies ( |Toft et ^|2012[ |2014) . Some of these objects appear 
to be the cores about which the most massive of today’s galaxies 
were built ^van Dokkum et al.|2014^ . Since evolved 2 ~ 2 galax¬ 
ies can serve as observational proxies for similar objects at higher 
redshifts, these interlopers are themselves promising candidates for 
future study. 

However, our primary interest for this work is high-redshift 
galaxies; thus, to exclude low-redshift objects we fit galaxy tem¬ 
plate spectral energy distributions (SEDs) to the observed fluxes to 
obtain photometric redshifts for our candidates. 

For this purpose we use BP zlj (Bayesian Photometric Red- 
shift: |BeMtez|[2000[ [Benitez et al.|20()4l [Coe et aL||2006| l, which 
matches an object’s SED to known galactic spectral types and pro¬ 
duces a probability distribution for the object’s redshift. 

However, BPZ’s default templates do not allow one to include 
extinction as a separate parameter, although they do empirically re¬ 
produce the observed photometry of the wide variety of galaxies 


^ We exclude detections within fifteen pixels of the edge of the frame as 
well as in the defect near the bottom of the detector (dubbed the ‘death star’ 
in |Dressel|2014j , in addition to overexposed regions, defined as any region 
with a surface brightness of less than 15 magnitudes arcsec“^. 

® http://www.stsci.edu/-dcoe/BPZ/ 
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Figure 2. Contours showing prior probability distribution for intrinsic ex¬ 
tinction Ay at a given redshift z, derived from data in|Bouwens et al.|p009| 
and jMagdis et al.|j2010) . 


observed in large-scale surveys (see, for example, [Rafelski et al.| 
|2015| >. We deemed it advisable to explicitly include a prior for in¬ 
trinsic exinction, given the incidence of galaxies with non-neglible 
dust attenuation (e.g., |Boone et al.|201 l[|Dey et al.|1999|>. We thus 
obtain dust-extinguished templates by applying the Calzetti et al.| 
( |2000^ extinction law to each default template in increments of 
AAv = 0.5 from Ay = 0 to 3. Not all extinction values oc¬ 
cur with equal probability, however, and the probability of extinc¬ 
tion evolves with redshift. We therefore require a Bayesian prior 
for the likelihood of a given extinction Av at a given redshift. To 
obtain such a prior we employ the data from fig. 6 of|B ouwens] 
|et al.| ( [200^ , who use observed UV continuum slopes /3 for a sam¬ 
ple of Lyman break galaxies to derive estimates of E{B — V) for 
redshifts from 2.5 to 6. |Bouwens et al.H2009t also estimate the ef¬ 
fective selection volume at each redshift, allowing calculation of 
the percentage of galaxies with a given intrinsic extinction (as a 
function of redshift). We extend these probabilities to the local uni¬ 
verse using [Magdis et ar]j2010^ ’s estimate that galaxies are 8 to 
10 times less obscured at z ~ 2 than they are now; and we ex¬ 
tend them to higher redshifts by assuming the same probabilities 
obtained at 2 : ~ 6. We then interpolate, smooth, and normalize the 
resulting distribution to obtain the prior shown in Fig.|^ One can 
obtain an analytic estimate for this distribution (Repp and Ebeling, 
in preparation), but for this work we simply employ the numerical 
probabilities plotted in the figure. 

In addition to the extinction prior, we also require some as¬ 
sumptions about the probability of observing various galactic spec¬ 
tral types at various redshifts. BPZ derives its default P{z\T) prior 
from the Hubble Ultra-Deep Field; however, since our area is so 
much larger than that of the HUDF, and since it - by design - in¬ 
cludes large-scale structure, one must adopt a broader prior that 
considers galaxies over a wider range of redshifts and masses than 
those encountered in the HUDF. In formulating an alternative to 
the default prior, we must balance the need to rule out inherently 
improbable results (such as observing an elliptical galaxy at 2 = 8) 
with our limited knowledge of galaxy evolution. 

We first rule out high-redshift ellipticals. We construct a 


Schechter luminosity function for elliptical galaxies using the fol¬ 
lowing parameters, derived by [Nakamura et al.| (|2003|( from the 
Sloan Digital Sky Survey (SDSS); M*(r*) = -21.52; a = 
-0.83; and (j>* = 1.61 x 10"^ Mpc"^ We then take the BPZ el¬ 
liptical galaxy template, redshift it, account for intergalactic atten¬ 
uation ( |Madau|1995^ , and convolve it with the F140W filter profile. 
Thus we derive an apparent magnitude for each combination of ab¬ 
solute magnitude and redshift. By combining this information with 
our F140W limiting magnitude and the elliptical galaxy luminosity 
function, we obtain a prior for observing elliptical galaxies which 
vanishes smoothly around 2 = 3.5. 

This constraint on elliptical galaxy visibility follows directly 
from the weakness of their UV emission combined with their em¬ 
pirically determined maximum luminosities. In recognition of the 
fact that the luminosity functions evolve with cosmic time (e.g., 
[Bouwens et al.|[2011| |2012^ , we assume flat priors for the other 
galactic spectral types. In particular, we have modified BPZ so that 
for each redshift, it reports the goodness-of-fit probability for the 
most likely galaxy template only, rather than summing the proba¬ 
bilities over all templates. (See also the discussion in Section HD 
concerning comparison of BPZ results with results from hyperz.) 
Thus we consider only how closely the observed SED fits a galaxy 
template at high redshift without taking into account the (unknown) 
likelihood of each template. 

To validate the utility of this modification, we next determined 
which procedure (original or modified BPZ) best reproduces the re¬ 
sults of CLASH using only our four passbands. Since CLASH uti¬ 
lizes 16 passbands for their photometric redshift determinations, 
their multiple bands in essence function as very low-resolution 
spectroscopy. Thus, reproducing their redshifts would enhance our 
confidence that our procedure yields reliable results. 

We began with 47 high-redshift galaxies from tabs. 5 and 6 
of |Bradley et al.|(2014| l. For each galaxy, [Bradley et al. [ report the 
photometric redshift estimate along with the 95 per cent (2cr) con¬ 
fidence interval. We obtained the magnitudes of each galaxy in our 
four passbands from the CLASH source catalog We then ran both 
the original implementation of BPZ and our modified version (us¬ 
ing the aforementioned prior on elliptical galaxies in both cases) 
on these 47 objects using only the four passbands considered in 
this project. From the BPZ results, we determined the most likely 
redshifts and 68 per cent (Icr) confidence intervals. Since we only 
use four passbands whereas CLASH uses 16, we expect their 2(j 
confidence intervals to be roughly comparable to our la intervals. 
Our limited number of passbands means that, in many cases, 95 per 
cent confidence intervals derived from our work would be so broad 
as to be almost useless. Thus for our own results we quote 68 per 
cent intervals and, as noted in Sectionj^ recommend spectroscopic 
follow-up for our most plausible candidates. 

The results appear in Fig.[^ We see that both versions of BPZ, 
when operating on only four passbands, give highest-likelihood 
redshifts roughly comparable to those obtained by CLASH. The 
original version of BPZ, as might be expected, produces a more 
symmetric scatter about the diagonal 24bands = 2 ;clash; it also 
results in fewer ‘catastropic outliers.’ However, the results of the 
modified version are more conservative in that they seldom produce 
redshifts significantly in excess of the CLASH results. The great¬ 
est excess redshift (compared with CLASH) is A 2 = 0.5 for the 
modified version, as opposed to A 2 =1.5 for the original version. 

In two other respects our approach is conservative: first, we 


https://archive.stsci.edu/prepds/clash/ 
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CLASH results (95% confidence) 


Figure 3. Comparison of photometric redshifts derived from CLASH with 
those derived from four passbands with BPZ (both the original version and 
as modified - see Section]^. Green dotted lines are Az = 1 away from 
the main diagonal Z 4 bands = -^^CLASH- The results of the original BPZ are 
more symmetric about the diagonal—and have fewer catastrophic outliers. 
However, the modified version produces more conservative results, in that 
it is less likely to produce photometric redshifts significantly in excess of 
the (presumably more nearly correct) CLASH estimates. 


apply our extinction prior to all spectral types of galaxies, thus 
assuming that dusty ellipticals are as likely as dusty starbursts or 
spirals. Second, we use the SDSS z = Q luminosity function for 
ellipticals to estimate the likelihood of observing such galaxies at 
higher redshifts, thus neglecting evolution. As a result of these two 
assumptions, our estimated likelihood for dusty z ~ 2 ellipticals 
is probably higher than that found in the actual universe. Given 
our limited number of passbands, our primary concern is to obtain 
conservative photometric redshifts; we are thus willing to accept 
results which might underestimate the true redshift. Hence we use 
the modified version of BPZ, with the understanding that we prob¬ 
ably underestimate to some degree the number of galaxies in each 
high-redshift bin. 

4.2 Eliminating Stars and Sub-stellar Objects 

Section 1^ outlines our two-fold approach to eliminating stars and 
sub-stellar objects from our list of dropouts. The first aspect of our 
approach is the morphological criterion displayed in Fig.[T] any ob¬ 
jects in the light green portions of this figure we accept as galaxies. 
However, objects in the blue portions of the figure have ambiguous 
morphology, and for these objects we use BPZ to determine how 
well their photometry matches what one would expect for M-stars 
and brown dwarfs. 


To do so, we obtained composite spectra of M-, L-, and T- 


dwarfs by stacking 25 spectra jBurgasser et aL||2004| |2006 

[Bur- 

gasser|2007a|b||Burgasser et al.|2008||2010||Chiu et al.|2006| 

Kirk- 

Patrick et al.|2010) 
Spectral Librariesl 

of these objects obtained from the SpeX Prism 


We then prepared simulated stellar observations by convolv¬ 
ing the model spectra with the HST filter profiles, scaling to mag¬ 
nitudes typical of our candidates and adding uncertainties typical 
of our candidates. Experimentation showed that, when applied to 
these simulated stars, BPZ with the extra templates yielded a prob¬ 
ability spike at 2 = 0 that was typically at least 30 per cent of the 
maximum height of the distribution. 

Thus, to impose our spectral criterion, we ran BPZ (with the 
extra templates) on each of our ambiguous candidates and elimi¬ 
nated those for which the probability density at z = 0 was at least 
30 per cent of the maximum probability density. The remaining ob¬ 
jects we retained as galaxies. Finally, we removed the stellar tem¬ 
plates from BPZ’s template library and re-ran it on these retained 
objects to obtain the probabilities reported in Table|^and shown in 
Figs|^and|^ 


5 DISCUSSION 


5.1 SED-fit Codes 


There seems to be no ‘standard’ code for SED-fitting. [Brad l ey et al.| 
l |2014[ use BPZ, whereas Zheng et al. 1 2014) use a combination 
of BPZ and iSEDfit. [Atek et al.H2oi4t use hyperz; |Bowler| 
jet al.| ) (2014| l use LePhare; and jOesch et al.| i 2013[ l use ZEBRA 
and EAZY. Others i jMcLure et al.|20lT||2013[ [E lis et al.|2013| l use 
proprietary code. 

Since there is no accepted best SED-fit code (see jHildebrandtj 
jet al.|20f0| for a review and evaluation of 17 photometric redshift 
codes) we thought it useful to compare the results of the two codes 
with which we are most f amiliar, namely BPZ ( modified as de¬ 
scribed above) and hyper j^) |Bolzonella et al. 2000 l BPZ marginal¬ 
izes over a carefully selected set of galaxy templates (a feature 
which our modification largely circumvents) and naturally accom¬ 
modates a prior distribution for those templates. It does not natu¬ 
rally handle various levels of intrinsic extinction, requiring us to 
handle extinction as described in Sectionj^ On the other hand, hy¬ 
perz fits not only redshift but also extinction and metallicity; how¬ 
ever, it does not seem to accommodate a prior on these variables. 

Thus we compare the results of two approaches: the first is a 
modified BPZ with an extinction prior and a limited galaxy tem¬ 
plate prior to disallow high-redshift ellipticals; the second is (un¬ 
modified) hyperz. The output appears in Fig. We note first the 
high incidence of objects (blue points) to which the codes assigns 
a probability distribution which is bimodal at a 68 per cent level 
(meaning that the 68 per cent confidence region is a topologically 
disconnected set). The majority of objects for which the 68 per cent 
BPZ and hyperz results differ (open circles) are these bimodal ob¬ 
jects with poorly constrained redshifts. Second, as in the compar¬ 
ison with the CLASH results, we see that our modified version of 


http: / /pono . ucsd. edu/ ~adam/br own dwarf s/ 
spexprism/library.html We then included these spectra in 
the BPZ template list and, for these three spectral templates only, imposed 
a delta-function prior limiting them to z = 0. 

® http://webast.ast.obs-mip.fr/hyperz/ 
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6 Repp et al. 



Figure 4. Photometric redshifts produced by BPZ (modified as described 
in Section]^ and by hyperz. Blue squares indicate probability distribu¬ 
tions from BPZ which are bimodal at the 68 per cent level (i.e., the 68 per 
cent confidence regions are topologically disconnected), with only the most 
likely redshift plotted for each candidate. Red circles indicate unimodal dis¬ 
tributions (connected 68 per cent confidence intervals). Filled markers indi¬ 
cate BPZ results which are consistent (to 68 per cent) with hyperz results; 
open markers indicate inconsistent results. Dotted lines show a distance of 
A 2 : = 1 from the main diagonal. For readability, we show error bars for 
consistent and unimodal distributions only. We use only the red data points 
(redshift probabilities unimodal at 68 per cent) for the remainder of our 
analysis. 


BPZ is conservative in that it tends (with a few exceptions) to as¬ 
sign a lower redshift to (unimodal) objects than does hyperz. Thus, 
consistent with our conservative approach, we shall henceforth ig¬ 
nore all objects to which BPZ assigns a bimodal distribution at the 
68 per cent level - i.e., we shall consider only the red points plotted 
in Fig. 12 


Table 3. Magnifications obtained for high-redshift candidates 


Object 

Photo- 

Ampli- 

Nominal 

iiQ 

1 M 1 

ficatioij^ 

Misocj^ 

MACSJ0140-0851 

7 

' ■'^-0.2 

1.26 ±0.00 

-22.8 ±0.3 

MACSJ0140-1028 

^•-*^-1.5 

7.04 ± 1.46 

-18.4 ±0.7 

MACSJO152-0651 

fi 7+0-4 
^•'-0.4 

10.03 ±0.26 

-18.6 ±0.2 

MACSJO152-0871 

6 

^•^-0.1 

7.25 ± 0.03 

-21.9 ±0.0 

MACSJ0712-0608 

7.0l°;» 

1.24 ±0.00 

-21.1 ± 0.2 

MACSJ0712-0699 


1.51 ± 0.01 

-23.1 ± 0.5 

MACSJ0947-0072 

fi ^ + 0-4 

1.17 ± 0.00 

-22.2 ± 0.1 

MACSJl 133-0922 

7 9+^’^ 

+ ^-1.0 

1.28 ±0.00 

-21.3 ±0.2 

MACSJ2135-0509 

6 9+° '^ 
O-y_0.4 

5.15 ±0.10 

-19.7 ± 0.2 

MACSJ2135-0763 

Q 1+0-3 
^•^-0.5 

1.41 ± 0.00 

-23.8 ±0.1 


^ SeeTablel 
b 
c 


El 


68 per cent confidence intervals. 


Derived from Lenstool ' Julio & Kneib 


2009 


models. 


Uncertainties propagated from photo-z uncertainties. 


( |2011^ note that UV magnitudes below —24 are physically un¬ 
likely, given that the vigorous star formation required for such lu¬ 
minosity would quickly fill the galaxy with dust. One of our objects 
(MACSJ2135-0763) approaches this limit of —24. One possible ex¬ 
planation is galaxy-galaxy leasing, given the nearby lower-redshift 
galaxy visible immediately to the left (in the third row of Fig.[^. 
Another possibility is that the photometric redshift, despite being 
constrained by BPZ at the 68 per cent level, is in error. Inspection of 
its redshift probability distribution (Fig.[^ shows a non-negligible 
probability of z ~ 2; nevertheless, the probability of z ^ 6.5 is 
around 85 per cent (see Table |^. This object would merit addi¬ 
tional investigation. 


5.2 Magnification 


Determining the rest-frame UV magnitude of dropouts requires 
knowledge of the degree to which gravitational leasing has mag¬ 
nified the object in question. Determination of the leasing magni¬ 
fication requires in turn knowledge of the mass distribution of the 
foreground cluster. 

Mass maps based on spectroscopy of strong-leasing features 
are available for a fraction of our target clusters, i.e. those for 
which such features have been spectroscopically confirmed. Mag¬ 
nification estimates were derived by means of Lenstool (|Jullo| 
|& Kneib]|2009| ) from the positions of strong-lensing features and 
a parametric model which includes the contribution from cluster- 
scale and galaxy-scale haloes (following the same approach as, 
e.g., [Richard et al.|2010[[Limousin et al.|2015^ . The available mod¬ 
els allowed us to calculate the magnification due to cluster gravity 
for the ten high-redshift candidates presented in Table[^ Note that 
the magnification errors include the statistical error from the para¬ 
metric model but no systematic uncertainties from different mod¬ 
elling approaches. Although the uncertainties listed in Tablej^also 
include errors propagated from redshift uncertainties, they should 
thus be regarded as lower limits. 

Some of these objects appear to be extremely luminous. Mag- 
nitudes below —22 are not unheard of|^ however, Bouwens et al. 


An example would be the object al a = 325.0753169°, (5 = 


6 RESULTS 

Given our extinction prior, our galaxy type prior, and our modified 
version of BPZ, we determine both a posterior probability distri¬ 
bution and a best-fitting SED for each of the dropout galaxies. We 
retain only the 50 objects which meet the following criteria. First, 
the BPZ-derived probability distribution must peak at 2 > 5.5. 
Second, the probability distribution must be unimodal at the 68 per 
cent level; by this we mean that the 68 per cent confidence region 
(shaded blue in Figs.|^and|^ must be a connected set, so that only 
one peak rises to 68 per cent significance. Table lists all 50 of 
these / 814 -dropout galaxies, their magnitudes (corrected for Milky 
Way extinction) in each of the four passbands; their photometric 
redshifts; and their probabilities of falling into redshift bins 2 ~ 7, 
2 ~ 8, and 2 ~ 9. 

Note that two of these objects (EMACSJ1057-2279 and 
EMACSJ1057-2476 in Table|^ appear to be multiply lensed im¬ 
ages of the same source. Both of these objects appear in the top 
row of Eig. 1^ the first is in the centre of the stamp; the other is 
barely visible at the right-hand edge of the stamp. We make no at¬ 
tempt to determine the number of other objects in our sample which 
might be multiply lensed, and in this work we do not account for 


—23.6772512° in table 4 of |Bradley et al.|(2014) , which yields arestframe 
Misoo around —23, although the authors note a small possibility that the 
unresolved source might be a star. 
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Object ID: EMACSJ1057-2279 
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Figure 5. Images (5 arcsec per side), redshift probability distributions, and best-fitting SEDs for selected dropout galaxies. For object IDs see Table]^ Blue 
shading denotes 68 per cent confidence (also noted above the probability plots). Red SEDs show the most probable fit; grey SEDs show the most probable 
low-redshift fit. The best-fitting template for each case appears above the SED plot. Green dotted lines show limiting magnitude in each filter. 
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Object ID: MACSJ1115-0632 
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Figure 6. Images (5 arcsec per side), redshift probability distributions, and best-fitting SEDs for selected dropout galaxies. For object IDs see Table]^ Blue 
shading denotes 68 per cent confidence (also noted above the probability plots). Red SEDs show the most probable fit; grey SEDs show the most probable 
low-redshift fit. The best-fitting template for each case appears above the SED plot. Green dotted lines show limiting magnitude in each filter. 


© 0000 RAS, MNRAS 000,pTfT2| 


































































































MACS High-Redshift Galaxies 9 


Table 4. Surveys plotted in Fig.^ 


Survey/ Depth*^ Solid Angle 

Reference Field(s) (5cr) (arcmin^) 


Atek et al. 

201 

4 

Abell-2744 (HFF) 

28.5 

4.7 

Bouwens et al. 

2011 

HUDF09, ERS 

28.5 

53 

Bouwens et al. 

2014 

CLASH 

27.4 

77 

Bradley et al. 

>014 

CLASH 

27.4 

82 

Fltis et at. 


ir 

HUDF 

29.5 

4.7 

Zheng et al 

2014 

Abell-2744 (HFF) 

28.5 

4.7 

This work 



MACS 

26.6 

135 


“Limiting magnitude in F140W filter, or inteipolated between F125W and 
F160W limiting magnitudes if F140W not utilized. If depths differ across 
fields in the same survey, we cite an average depth weighted by the solid 
angle of the fields involved. 

^In image plane, not source plane. 


w 1 
C 10“' 



Atek 

Bouwensll 

Bouwensl4 

Bradley 

Ellis 

Zheng 

This work 


27.5 28.0 28.5 29.0 

Limiting F140W magnitude (5a) 


29.5 


30.0 


this source of systematic error in our sums of objects in each red- 
shift bin. 

The FI low photometry of four of these objects is defective, 
as noted in Section]^ The objects appear at the end of Tablein 
addition, one of them appears in the fourth row of Fig. and an¬ 
other appears in the fourth row of Fig. For these objects, BPZ 
yields a relatively flat probability distribution at higher redshifts; in 
such cases, we weight the probability of galaxies’ placement into 
redshift bins with the summed probabilities for the other dropout 
galaxies - in essence using the other galaxies’ redshifts as a prior 
for those with defective FI lOW photometry. Figs.[^and|^show the 
images, probability distributions, and SEDs for some of the galax¬ 
ies with the greatest likelihood of lying at a high redshift. 

By summing the probabilities in Table we obtain the total 
number of galaxies detected in each bin (also reported in Table[^, 
which round to 19 galaxies at 2 ; ~ 7, 6 at 2 ; ~ 8 , and 3 at 2 ~ 9. 
We thus estimate a total number of 27 high-redshift galaxies, which 
is one less than the sum of the above numbers due to rounding. 
(See the table for 68 per cent Poisson confidence intervals.) These 
sums include one galaxy (MACSJ1319-0075) for which the BPZ 
confidence interval 2 = Q.SIq ® is inconsistent with the hyperz 
confidence interval 2 = 5. 7^3 5 . This object corresponds to the 
single open red circle in Fig. and appears with a footnote in the 
table. 

Finally, in order to evaluate the effectiveness of our approach, 
we compare the number of high-redshift galaxies estimated in this 
survey with those detected by others. One would expect a positive 
correlation between the number of galaxies and the limiting mag¬ 
nitude of the survey. Using a thousand 0.4"-diameter apertures in 
each of our clusters, we determine that the 5a limiting magnitude 
for our survey in F140W is 26.6. Normalizing the number of de¬ 
tected high-redshift objects to the solid angle coverage of various 
searches, and plotting the results against limiting magnitudes, we 
obtain Fig.|7] (See also Table|^for the fields involved.) 

Note that this figure combines surveys with disparate param¬ 
eters (e.g., various detection bands, differing uniformity of depth, 
lensed/unlensed, etc.); thus one ought not to attempt to derive any 
sort of cumulative luminosity function from it. However, what 
Fig. 0 does show is that our survey returns a number of galax¬ 
ies commensurate with our low limiting magnitude. It is the large 
solid angle of our survey that permits discovery of a non-negligible 
number of high-redshift candidates, underlining that shallow, wide 
surveys efficiently probe the bright end of the high-redshift UV lu¬ 
minosity function. 


Figure 7. Compaiison of the number of objects found in this work to that 
detected by other surveys, as a function of limiting magnitude. Filled sym¬ 
bols denote lensed surveys, and unfilled symbols denote unlensed surveys. 
(See Tablej^for details.) Blue markers indicate objects at 2 ~ 7; green, at 
2 ~ 8; and magenta, at 2 ~ 9. Markers for 2 ~ 8 and 2 ~ 9 are slightly 
offset horizontally for clarity. Error bars represent Poisson confidence in¬ 
tervals. Note that for lensed surveys the vertical axis lists the number of 
objects per unit solid angle in the image plane, not the source plane. 


7 CONCLUSIONS 

We draw the following conclusions from our survey. First, broad 
and relatively shallow surveys like ours can effectively detect high- 
redshift galaxy candidates. Devoting only about an hour of Hubble 
time to each of 29 clusters, we were able to identify 50 objects 
which have a probability distribution which peaks at high redshift 
and is unimodal at the 68 per cent level. Of these objects, we es¬ 
timate that 27 lie at redshift z > 7. Since our survey is relatively 
shallow, we are therefore sampling the bright end of the UV galaxy 
luminosity function at such redshifts. Our currently available mass 
maps suggest the existence of several extremely luminous objects 
with a photometric redshift of 2 ~ 8-9, placing them only about 
600 Myr after decoupling. Such objects are ideal targets for fur¬ 
ther study, if for no other reason than to confirm or falsify the high 
photometric redshifts. 

This further investigation is important for at least three rea¬ 
sons. First, despite our rejection of objects with significantly bi- 
modal distributions, inspection of Table shows that almost all of 
the remaining objects still have non-negligible lower-redshift prob¬ 
abilities. Indeed, the summed likelihoods show that we can expect 
only 27 of the 50 listed objects to actually lie at 2 > 7. Second, 
we derived our extinction prior from Lyman-break galaxy data, 
whereas low-redshift F814W dropouts represent a different popula¬ 
tion (4000-A break galaxies and possibly some high-luminosity in¬ 
frared galaxies). This population mismatch at lower redshifts could 
potentially introduce a systematic bias into our results. Third, our 
analysis of BPZ and hyperz shows that the two codes can occa¬ 
sionally yield wildly divergent results. Thus the choice of a specific 
SED-fit routine can introduce additional - and probably significant 
- systematic bias. We are unable to quantify these biases, and the 
resulting intrinsic uncertainty is not captured in the statistical error 
bars we report. 

We note, however, that the low-redshift interlopers are scien¬ 
tifically interesting in their own right. Over 70 of our 124 F814W- 
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dropout galaxies do not appear in Table due to bimodalities in 
their probability distributions; thus, besides high-redshift candi¬ 
dates like those we report in this work, surveys such as ours can 
identify a large number of lower-redshift objects conducive to the 
study of galaxy evolution. 

Finally, we conclude that even as few as four hroad-band fil¬ 
ters suffice to isolate promising objects. However, our results also 
reinforce the importance of SED fitting in obtaining reliable pho¬ 
tometric redshifts: dependence on a simple dropout criterion alone 
would have doubled our catalog size by introducing objects with a 
significant likelihood of lying at low redshift. SED-fitting is espe¬ 
cially important for surveys which utilize a small number of pass- 
bands; a greater number of passbands allows the use of more com¬ 
plex colour criteria (e.g., [Castellano et al.|20ldl l which in practice 
function as a coarse SED fit. In particular, the F606W band was 
necessary to eliminate high equivalent-width emission-line galax¬ 
ies from our sample. The small number of passbands did lead to 
a large number of objects with significantly bimodal photometric 
redshifts, which we eliminated from consideration; it was the large 
solid angle of the survey which permitted isolation of a significant 
number of more securely identified objects. 

Spectroscopic investigations of some of our best candidates 
promise to better constrain these objects’ redshifts, as do deeper 
space-based observations with more filters, and/or slitless near in¬ 
frared spectroscopy with WEC3 grisms. In addition, the eMACS 
SNAPs program proposes to survey 50 massive clusters with z > 
0.5. Application of our method to this and other future programs 
has the potential to isolate yet more promising candidates. Thus we 
confidently expect that the study of high-redshift galaxies will con¬ 
tinue to shed light on the processes occurring in the early universe. 
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Table 5. Dropout galaxies; all confidence intervals are 68 per cent. 


ID 


“J2000 


<5j2000 


”^F606W 


«tF814W 


t«F110W 


t«F140W 


Ps® Pg® 


EMACSJ1057-226I 

10*'57™32.02= 57°59'30" 

> 28.8 

27.48 ± 0.69 

26.04 ±0.13 

26.17 ± 

0.20 


0.19 

0.01 

0.00 

EMACSJ1057-227^ 

10'’57™27.80= 

57° 59'7" 

> 28.8 

> 28.4 

25.15 ± 0.06 

24.86 ± 

0.06 

6-9la3 

0.79 

0.14 

0.00 

EMACSJ 1057-2476^ 

10‘'57™27.48= 

57° 59'7" 

> 28.8 

> 28.4 

24.55 ± 0.03 

24.15 ± 

0.03 

6.9lg;f 

0.74 

0.25 

0.00 

MACSJ0140-0851 

1*'39“58.63= 

-5°56'12" 

> 28.8 

> 28.4 

25.13 ±0.05 

24.25 ± 

0.03 

7.8l°;l 

0.19 

0.60 

0.17 

MACSJO140-1028 

1*M0™3.70= 

-5°55'15" 

> 28.8 

28.04 ±0.97 

26.38 ±0.19 

26.25 ± 

0.24 

'^•-‘-- 1.8 

0.19 

0.09 

0.01 

MACSJ0152-0477 

1*'52“33.16= 

-28°54'42" 

> 28.8 

> 28.4 

26.62 ± 0.16 

26.35 ± 

0.17 


0.32 

0.17 

0.01 

MACSJO 152-0651 

1*'52“34.36= 

-28°54'27" 

> 28.8 

> 28.4 

25.89 ±0.11 

25.88 ± 

0.16 

6 - 8 lg:i 

0.67 

0.07 

0.00 

MACSJ0152-0871 

1*'52“36.05= 

-28°54'35" 

> 28.8 

> 28.4 

22.76 ±0.02 

22.79 ± 

0.03 

6-9lg:J 

1.00 

0.00 

0.00 

MACSJ0257-0913 

2h57“41.18= 

- 22 ° 10 ' 10 " 

> 28.8 

> 28.4 

26.23 ±0.15 

25.89 ± 

0.14 

6-9lo:9 

0.37 

0.23 

0.01 

MACSJ0712-0608 

7*'12™13.17= 59°32'52" 

> 28.8 

> 28.4 

25.99 ±0.10 

25.58 ± 

0.10 


0.43 

0.25 

0.00 

MACSJ0947-0072 

9*'47“17.86= 76°24'21" 

> 28.8 

27.87 ± 1.37 

24.62 ± 0.05 

24.56 ± 

0.07 

6 d 
+ 1 

CO 

0.53 

0.00 

0.00 

MACSJl 115-0329 

11*M5™15.94= 

53°19'5" 

28.28 ± 0.84 

> 28.4 

26.35 ± 0.18 

26.34 ± 

0.25 

6 . 811:1 

0.39 

0.13 

0.01 

MACSJl 115-0632 

ll‘'15™17.45= 53°20'27" 

28.65 ± 1.60 

> 28.4 

25.23 ±0.08 

24.67 ± 

0.07 

7 O+0.7 
'•'^-0.7 

0.48 

0.39 

0.01 

MACSJl 124-0811 

11*’24™28.77= 43°50'41" 

> 28.8 

> 28.4 

26.38 ±0.17 

25.89 ± 

0.15 

y 1 +1-3 
' •-‘--1.4 

0.28 

0.25 

0.04 

MACSJl 133-0922 

Ilh33"i8.33= 

50°8'27" 

28.63 ± 1.39 

> 28.4 

25.98 ±0.13 

25.47 ± 

0.12 

7-211:1 

0.35 

0.29 

0.02 

MACSJl 621-0860 

16*'21™23.01= 38°11'13" 

> 28.8 

> 28.4 

26.51 ± 0.13 

26.22 ± 

0.14 

6 . 811:1 

0.36 

0.17 

0.00 

MACSJl 652-0135 

16*’52™26.52= 55°34'43" 

> 28.8 

> 28.4 

26.81 ± 0.19 

26.52 ± 

0.20 

6 9^^'^ 

’^•^-1.9 

0.26 

0.17 

0.02 

MACSJ2051-0806 

20*'51™13.80= 

2°16'48" 

> 28.8 

> 28.4 

26.32 ± 0.16 

26.11 ± 

0.18 

'^•^- 1.0 

0.39 

0.18 

0.01 

MACSJ2135-0509 

21*’35"'10.90= 

-1°3'12" 

> 28.8 

> 28.4 

25.76 ±0.14 

25.52 ± 

0.16 

q+0.9 

0.53 

0.23 

0.00 

MACSJ2135-0763 

2l''35“8.24= 

-1°2'41" 

> 28.8 

27.80 ± 1.11 

24.59 ±0.04 

23.36 ± 

0.02 

Q 1 +0-3 

^•-‘^-0.5 

0.00 

0.07 

0.78 

MACSJ2135-1078 

2lh35™14.82= 

-1°2'25" 

28.72 ± 2.05 

> 28.4 

26.15 ± 0.18 

25.69 ± 

0.16 

7 9 + 1-2 

0.31 

0.26 

0.05 

SMACSJ0600-0180 

6*'0“9.21= 

-43°53'6" 

> 28.8 

28.04 ± 1.19 

25.48 ±0.11 

24.81 ± 

0.08 


0.27 

0.26 

0.03 

SMACSJ0600-0427 

6*'0“9.12= 

-43°53'44" 

> 28.8 

28.33 ± 1.44 

25.01 ± 0.07 

24.06 ± 

0.04 

0 ri + 0.9 
^■*-'-0.5 

0.11 

0.48 

0.25 

SMACSJ2031-0768 

20*’31™45.98= 

-40°37'31" 

> 28.8 

> 28.4 

26.38 ±0.13 

25.87 ± 

0.11 

7-211:1 

0.30 

0.25 

0.02 

SMACSJ2131-0444 

2l''31“6.95= 

-40°18'57" 

> 28.8 

> 28.4 

24.38 ± 0.04 

24.33 ± 

0.06 

6.9l°:l 

0.98 

0.00 

0.00 

SMACSJ2131-0516 

21>'31“6.47= 

-40°18'51" 

> 28.8 

> 28.4 

25.07 ±0.08 

24.83 ± 

0.09 

6.9l°:® 

0.78 

0.16 

0.00 

SMACSJ2131-0567 

2lh3ini6.14= 

-40°18'23" 

27.61 ± 0.91 

> 28.4 

25.24 ±0.09 

24.45 ± 

0.06 

7 K + 1-1 
‘ •^-0.4 

0.35 

0.46 

0.11 

Summed probabilities for objects with reliable photometry in 

all four bands: 






11.6 

5.6 

1.6 

Dropouts with defective FllOW photometry 

0 










MACSJ0712-0414 

7*'12“25.64= : 

M°31'55" 

> 28.8 

> 28.4 

— 

22.72 ± 

0.02 

Q 9 + 1.8 

0.21 

0.31 

0.31 

MACSJ0712-0699 

7*'12™29.23= 59°32'59" 

> 28.8 

27.20 ± 1.55 

— 

23.69 ± 

0.05 

^ [-+ 2.0 

0.18 

0.17 

0.17 

SMACSJ0549-0900 

5*'49“16.25= 

-62°5'15" 

> 28.8 

> 28.4 

— 

24.35 ± 

0.05 

8 -oll:| 

0.23 

0.24 

0.24 

SMACSJ0549-1009 

5h49m;^4g7s 

-62°5'46" 

> 28.8 

> 28.4 

— 

23.68 ± 

0.04 

7-911:2 

0.23 

0.28 

0.28 

SMACSJ0549-9147 

5*'49“11.29= 

-62°5'7" 

28.56 ±2.22 

> 28.4 

— 

24.62 ± 

0.07 

7-911:1 

0.22 

0.23 

0.23 


Summed probabilities for candidates with defective FI lOW photometry p| 
Net summed probability in each redshift bin: 

Net detections reported in each redshift bin: 


0.7 

12.2 

12 


+4 

3 


0.4 

5.9 

6+3 

”-2 


0.1 

1.7 

2 


+ 2 
-1 


“ 68 per cent confidence intervals. 

^ P 7 , Pg, and P 9 denote the probabilities that the object falls within redshift bins [6.5, 7.5), [7.5, 8.5), or [8.5, 9.5), respectively. 

^ Part of multiply-lensed system; see Sectionj^ 

Reported probabilities for dropouts with defective FllOW photometry are those derived directly from the BPZ fit, before application of a prior based on the 
redshifts of the other dropouts (see text). 

® Summed probabilities reflect the prior described in the text and referenced in the immediately preceding note. 
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